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We analyze the group-theoretical ramifications of the Nambu-Goldstone [NG] 
(^T) ' theorem in the self-consistent relativistic variational Gaussian wave functional ap- 

proximation to spinless field theories. In an illustrative example we show how the 
Nambu-Goldstone theorem would work in the O(N) symmetric tp 4 scalar field the- 
ory, if the residual symmetry of the vacuum were lesser than O(N-l), e.g. if the 
vacuum were 0(N-2), or 0(N-3),... symmetric. [This does not imply that any of 
the "lesser" vacua is actually the absolute energy minimum: stability analysis has 
not been done.] The requisite number of NG bosons would be (2N - 3), or (3N - 
6), ... respectively, which may exceed N, the number of elementary fields in the La- 
grangian. We show how the requisite new NG bosons would appear even in channels 
that do not carry the same quantum numbers as one of N "elementary particles" 
(scalar field quanta, or Castillejo-Dalitz-Dyson [CDD] poles) in the Lagrangian, i.e. 
in those "flavour" channels that have no CDD poles. The corresponding Nambu- 
Goldstone bosons are composites (bound states) of pairs of massive elementary 
(CDD) scalar fields excitations. As a nontrivial example of this method we apply 
it to the physically more interesting 't Hooft a model (an extended Nf = 2 bosonic 
linear a model with four scalar and four pseudoscalar fields), with spontaneously 
and explicitly broken chiral 0(4) x 0(2) ~ ST/r(2) x SU\,(2) x ?7a(1) symmetry. 



1 Introduction 

The proof of the Nambu-Goldstone [NG] theorem [1, 2, 3, 4, 5] in the Gaussian wave 
functional approximation [6, 7, 8, 9, 10, 11, 12, 13] used to be an open problem for 
over 30 years, see Refs. [14, 15, 16, 17]. The first proof was given in the 0(2) 
symmetric cj) theory [17], and then straightforwardly extended to 0(4) in [18]; the 
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crucial assumption was that the ground state be O(N-l) symmetric, i.e. that only 
one (scalar) field develops a vacuum expectation value [v.e.v.]. 1 By the standard 
NG boson counting methods [3, 4], for every spontaneously broken symmetry Lie 
group generator there is one NG boson. As the Lie algebra O(N) has N(N — l)/2 
generators, for a ground state (vacuum) with an O(N-l) residual symmetry the 
number of NG bosons ought to be N(N - l)/2 - (N - 1)(N - 2)/2 = N - 1. That 
is exactly the number of available fields in the Lagrangian. What happens when 
the residual symmetry the ground state is "smaller" than O(N-l) and there ought 
to be more than N-l NG bosons? 

In this paper we shall extend our proof of the NG theorem in the Gaussian 
approximation [17, 18]to the O(N) symmetric <ft 4 theory when the symmetry of the 
ground state is dynamically broken to some (proper) subgroup of O(N-l), in this 
specific case to one of the following symmetries: 0(N — 2) x 0(2), 0(N — 3) x 
0(3), ... , 0(N/2) x 0(N/2) for N even, or 0((N - l)/2) x 0((N + l)/2) for N 
odd. As we shall show, the residual symmetry pattern is dictated by the absolute 
minimum of the 0(N) <fi 4 model's ground state energy, which in turn depends on 
the dynamics (gap equation) in the Gaussian approximation, as well as the free 
parameters. Absolute minimization of the vacuum energy will not been done in 
this paper, only a search for the extremal and/or saddle points. In a different 
model or approximation, or for different values of the free parameters the residual 
symmetry might be different. In any case, with such an asymmetric ground state 
the number of NG bosons must exceed N-l, the largest number of "elementary" 
scalar fields available in the Lagrangian (at least one field must develop the vacuum 
expectation value and thus cannot create or destroy single NG bosons). Nevertheless 
the canonical number of massless spinless excitations appears in the spectrum [4]. 
Our proof ought to leave no doubt as to the composite nature of the NG bosons in 
the Gaussian approximation. 

This paper falls into five Sections and two Appendices. Firstly in Sect. 2 we 
define the 0(N) </> 4 model and the Gaussian approximation. In Sect. 3 we show 
exactly how the requisite number of massless NG (bosonic bound) states appear and 
that all the "broken symmetry" Nother currents remain conserved, as the symmetry 
of the vacuum is reduced. NG bosons invariably appear in the Gaussian approx- 
imation to the 0(N) 4> 4 model in those channels that also contain CDD poles, so 
in Sect. 4 we extend our proof to the 't Hooft <f> A model which does not have this 
property. There we show how various bound NG states appear, or disappear as the 
symmetries of the ground state and/or the Lagrangian change. Finally, in Sect. 5 
we draw our conclusions and set them in a wider context. In the Appendices we 
present some technical details omitted in the main body of text. 

2 Preliminaries 

2.1 The O(N) symmetric scalar 4 model 

At first we confine ourselves to the 0(N) symmetric scalar 4 theory for the sake 
of simplicity. All scalar field theories with other spontaneously broken internal 

1 That assumption is justified when the remaining N-l fields are pseudoscalars, as in the Gell-Mann- 
Levy model, otherwise CP symmetry is (spontaneously) broken. 
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symmetries can be reduced to some subgroup of O(N). Of course, in such cases 
there will be interaction terms other than the simple 4> 4 one shown below. The 
Lagrangian density of this theory is 

c = lid^f-vitf), (1) 

where 

4> = (^0, 4>i, <h, <t>N-i) = 0, ?r), 

is a column vector and 

We assume here that Ao and Hq are not only positive, but such that spontaneous 
symmetry breakdown (SSB) occurs in the Gaussian approximation [GA] introduced 
below. 



2.2 The Gaussian variational method 
2.2.1 The Gaussian ground state ("vacuum") 

The Rayleigh-Ritz variational approximation to quantum field theories is based on 
the ("elliptical") Gaussian ground state (vacuum) functional Ansatz, Refs. [9, 10, 
11,13], 

= AAexp (-i- J dx| dy [&(x) - <&(x)>] G~\x,y) fo(y) - <&(y))>]) , (2) 

where N is the normalization constant, and one sums all repeated (roman lettered) 
indices from to N-l. (<pi(x)) is the vacuum expectation value (v.e.v.) of the 
z-th spinless field which henceforth we will assume to be translationally invariant 

<&(x)> = <<M0)> = (</>i) and 

G -U v) ~-S [ — 1 e 'k-(x-y) 

^(x,y)- 2 ^y (2vr)3y^^ e 



tir 



We have explicitly kept h (while setting the velocity of light c = 1) to keep track of 
quantum corrections and count the number of "loops" in our calculation. Then the 
"vacuum" (ground state) energy density becomes 



S(rrH,(<f>i)) = -^^ 2 + T [ (0)2 1 



2 

N-l r 



+ nj2 



i=0 



him) - 2 (Mo + m i ) I o( m i 



, (- N-l N-l 

+ 7 6h S (^) 2/ oK) + 2^ <^> 2 /oK 
I i=o «^i=o 

JV-l JV-1 

i=0 «^j=o J 



(3) 
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where 



t < ^ 1 f dk 1 . f cPk 1 _ , , fA . 

Ioim) = 2 J TPT^f = v wr i^^TTe = G ^ x) (4) 

h{rm) = \ J J^yJ k2 + m i =~lj 7^4 lo S ( fc2 " m ? + ie ) + const - ■ ( 5 ) 

We may identify KI\{rni) with the familiar "zero-point" energy density of a free 
scalar field of mass rrii . 

The divergent integrals Io i{mi) are understood to be regularized via an UV 
momentum cutoff A. Thus we have introduced a new free parameter into the cal- 
culation. This was bound to happen in one form or another, since even in the 
renormalized perturbation theory one must introduce a new dimensional quantity 
(the "renormalization scale/point") at the one loop level. We treat this model as 
an effective theory and thus keep the cutoff without renormalization 2 . 



2.2.2 The vacuum energy minimization equations 



We vary the energy density with respect to the field vacuum expectation values [fa) 
and the "dressed" masses rrii . The extremization condition with respect to the field 
vacuum expectation values (v.e.v.) reads: 



d£(mj, (4>j)) 



= 0; je0,l,...,7V-l ; 



(6) 



or explicitly 



N-l 



-vl + A ( <0> 2 + 3R7oK) + h ]T J ( 



rrii 



(7) 



The second set of energy extremization equations reads 
/ d£(rrii, {<pi})\ 



dm-i 



= 0, j€0,l,...,JV-l, 



(8) 



or 



2 i 2 

rrij + fiQ 



N-l N-l 

A ( 3<$> + h + 3WoK) + ?i E ^o(m fc ) 



(9) 



Eqs. (7), (9) can be identified with the (truncated) Schwinger-Dyson (SD) equation 
[13, 19] for the one- and two-point Green functions, see Ref. [17, 18]. The solutions 
to these equations plus the additional minimization requirements (positive-definite 
second derivatives matrix, i.e., the positive definiteness of its principal minors, which 
will not be discussed here) determine the symmetry of the ground state (vacuum). 



2 There are several renormalization schemes for the Gaussian approximation, but they show signs of 
instability and ultimately seem to lead to "triviality" [12]. 
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2.2.3 The vacuum symmetry 



Now one ordinarily assumes 3 that only one of the scalar fields <j>i [by convention 
the i = one] develops a nonzero v.e.v., i.e. (<fo) / and then one proceeds with 
the proof of the NG theorem. 

In the conventional case the scalar field masses are tuq = M; m\ = m-i = m% = 
... = n, i.e., the (N - 1) fields 4>i;i G 1,2,..., N — 1, of mass = /i, form an (N 
- l)-plet and the residual symmetry of the vacuum is O(N-l). That, however, is 
not the only logical possibility: one may assume that more than one field develops 
v.e.v.. This, of course, means that the (residual) symmetry of the vacuum is "lesser" 
than the one in the ("canonical") case with only one v.e.v. 4 . 

The decision which of these possibilities actually takes place can be made on 
the basis of comparing their respective ground state energies. In the case with an 
O(N-l) symmetric ground state, the vacuum energy, or effective potential is 



Similar expressions for ground states with lesser symmetries than O(N-l) can be 
derived by applying the corresponding gap equations (7), (9) to Eq. (3). The ques- 
tion of the absolute energy minimum will not be pursued in this paper, however. 
For the purpose of argument we shall assume that alternative ground states exist 
and are stable, i.e., energetically favourable to the standard one. 

If two (or more) scalar fields' v.e.v. are simultaneously nonzero, e.g. if (4>o) = 



Vl /Oand ( < f, j )=v 2 ^0;j € l,2,..,iV-l (</>■) =0;i € 1, 2, ...J - l,j - 2, .., TV - 1, 



we may change the scalar field labels such that the fields with nonzero v.e.v.s are 
labeled successively 0,l,...,k without loss of generality. In such a case the masses ra; 
may also take on different values. 

For simplicity's sake we assume that only two fields have nonvanishing vacuum 
expectation values (v.e.v.) 5 , 4>o,(f>i / (we label their masses mo = M,m\ = fii). 

3 The validity of this assumption, of course, depends on the values of the bare parameters in the 
Lagrangian and the cutoff A, or the values of the renormalized parameters, if one insists on renor- 
malization. One may also have an unbroken symmetry: with all the fields having a zero v.e.v., 
((pi) = 0,i £ 0,1, N — 1, the gap equations (7), (9) lead to all the masses rrii being equal. In other 
words the N fields form an O(N) multiplct, so we may say that the symmetry of the vacuum is O(N), i.e. 
not broken. 

4 In the Born approximation such a vacuum might be reducible to the "one v.e.v." vacuum by means 
of an O(N) transformation (if it lies in the same "orbit" of the symmetry group [20]), but in the Gaussian 
approximation dynamical symmetry breaking may lead to an irreducibly different ground state. 

5 Extension to three, or more v.e.v.s follows by straightforward analogy. 
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+ 



(10) 
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The remaining (N - 2) fields fai G 2, ...,(N — 1), with mass rrii = fi2, form an 
(N - 2)-plet. In that case, it is clear that the residual symmetry of the ground 
state is (at least) 0(N - 2), i.e., that the O(N) symmetry has been dynamically 
(spontaneously) broken to (at least) 0(N - 2). We say "at least" because, when the 
masses tuq = M = m\ = \x\ are equal there is additional 0(2) symmetry of the 
vacuum. Similar comments are valid in the case when more than two fields develop 
v.e.v.s. Of course, the residual symmetry of the vacuum determines the number of 
the NG bosons. Next we turn to the case of two fields with v.e.v. as worked out in 
Appendix A. We have shown in Appendix A that there are only two distinct mass 
solutions M = fi±, ^ [12 to Eqs. (55-58). Thus the residual symmetry of the vacuum 
is 0(N — 2) x 0(2), and the corresponding number of NG bosons must be (2N - 4). 
In the following we shall see how all these NG bosons come about. 

3 The Nambu-Goldstone theorem 
3.1 The two-body equation 

In Ref. [17, 18] it was shown that in the Gaussian approximation to the 0(2) and 
0(4) (f> 4 model, the Nambu-Goldstone particles appear in the two-particle spectra, 
i.e. that they are massless bound states of two different massive elementary ex- 
citations with an admixture of the (massive elementary) one-body state with the 
same quantum numbers (the CDD pole). This admixture of the one-body state is 
crucial for the masslessness of the NG state [it also proved to be a source of con- 
fusion]. As there are at most (N - 1) such elementary particles/CDD poles 6 , it 
appears that there can be at most (N - 1) NG bosons in this theory. There are 
N(N — l)/2 (distinct) pairs of elementary particles in this model, however, and an 
equal number of distinct (potentially bound) two-body states. Thus there can be 
at most N(N — l)/2 NG bosons, precisely the maximum number allowed by the 
0(N) Lie group generator counting. Next we must show that the number of "broken 
symmetry generators" corresponds precisely to the number of massless bound states 
and that the corresponding Nother currents are conserved. 

The two-body equation of motion in the Gaussian approximation is equivalent 
to the four-point Schwinger-Dyson (or Bethe-Salpeter) equation [for proof of this 
equivalence, see Ref. [9, 21]]. All the NG channels obey a generic four-point SD 
equation that reads 

(s) (11) 

where D(s) is the four-point Green function (N x N) matrix (scattering amplitude), 
n(s) is the polarization function matrix and V(s) is the potential matrix; ij denote 
the 0(N) indices of the two constituents, with the generic solution 

D( S )=V( S )[l-V( S )n( S r 1 (12) 

where s = (p± + p 2 ) 2 = P 2 is the center-of-mass (CM) energy. These matrices may 
reduce to a direct sum of submatrices depending on the residual symmetry of the 

6 0ne field adopts a v.e.v. so it cannot create or destroy single NG bosons. 
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system. Such effective two-body propagators can aiso be written in the foilowing 
form (see Ref. [22]) 



D a (s) 



s -mt: 



(13) 



where m a is the effective mass in channel a. The difference between various O(N) 
"flavour" sectors appears in the polarization functions Il(s) and potentials V(s): 
more specifically in (N - 2) channels [(i = 0,j C 2,...N — 1) = M//2] the said 
matrices are diagonal and have the form 



d 4 k 



1 



(2vr) 4 [k 2 — M 2 + ie] [{k - P) 2 -f4 + ie] 



Vmm ( s ) = 2A 



1 + 



2Aq(0q)- 
s-nl 



= 2A 


M 2 


1 + 2 




8-f4 



(14) 
(15) 



and another (N - 2) channels [(i = 1, j C 2, ...iV — 1) = fiifi2 have the form 

d A k 1 



/" d 1 



2A n 



1 + 



2A 



(2vr) 4 [k 2 -fi 2 + ie] [(k - Pf - fi 2 + ie] 

2" 



= 2A 









(16) 
(17) 



The poles 1 . 2 in Eqs. (15), (17) are called the Castillejo-Dalitz-Dyson [CDD] 

s /i 12 

poles and correspond to one-particle states in the theory. Now due to Eqs. (55-58) 
these two channels happen to be equivalent in this case, but that degeneracy is 
accidental (for a different example see Sect. 4). 



3.2 Massless (Nambu-Goldstone) two-body states 

With an 0(N — 2) x 0(2) residual symmetry there ought to be (2N - 4) = (N - 
2) + (N - 2) NG bosons. This number of NG bosons quickly exceeds the number 
of available fields in the Lagrangian: 2(N - 2) > (N - 1) for N > 3. Thus, there 
aren't enough scalar fields to provide for all the NG bosons in case of nonstandard 
symmetry breaking. Does this mean that the NG theorem breaks down in such a 
case? 

As in Ref. [17], we shall prove that at zero CM energy P = 0, the matrix 
[1 — ^2(0)^2(0)] vanishes. We use Eq. (49) to write 



1 - 



A*2 



V M ^ 2 (0) = 2A 

n MM2 (0) - hi — —2 — -2 — 



v m 2 - ^2 

then use Eq. (33) to obtain the final result 

VA^ 2 (o)n MAt2 (o) = 1 



(18) 



(19) 
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Thus the inverse propagator Eq. (12) evaluated at zero momentum vanishes: 



D 



(0) 



[i - yM M2 (o)n M ^(o)]VM^ 2 (o)- 1 = o, 



(20) 



which is equivalent to the result rnj^ = 0. Similarly for the remaining (N - 2) 
channels, q.e.d. 

Note that the proof in no way depended on the fact that /xi = M. In other 
words, if the gap Eqs. (55-58) allowed a solution such that fx± ^ M and {(f>\) = 0, 
the NG theorem would still hold. But, then the CDD poles would decouple in the 
latter (N - 2) channels, i.e., there would be no CDD poles in these channels. This 
shows that, at least in principle, (N - 2) NG bosons could be pure bound states with 
no single particle admixtures. In Sect. 4 we shall give an example of a model with 
a gap equation that allows such solutions. 



3.3 Conservation of Nother currents 

There are a C 2, ...,N — 1 dynamically broken O(N) symmetry Nother current 
matrix elements corresponding to 

JU(p>p) = (r(p)\Mo)\Mp)) 

- TjS*{q)D Mlta {q) . 



(21) 



where T^(q) is 



£t5 



ta) - tf 



„2 



d A k 

(2vr) 4 

2 



(2k + q)„ + q„ 



M 2 



q 2 -iA 



[k 2 - M 2 ] [(k + q) 2 



^- (vA^ 2 (o)n MM2 (o) - v M ^(q 2 )n Ml , 2 (q 2 )) 



(22) 



Inserting the vertex T^ 2 (g), Eq. (22) together with the two-body propagator 
D M ^ 2 {q 2 ), Eq. (12) into Eq. (21) one finds 

J a ^{p\p) = ip' +P)» + q» ( ^2^2 j 



(23) 




where q u = (p — p) u . This current is manifestly devoid of a pole at q 2 = fi 2 . 
The composite state plays precisely the role of the Nambu-Goldstone boson in the 
conservation of the dynamically broken O(N) symmetry Nother currents [1], i.e., in 
the basic O(N) symmetry Ward-Takahashi identity, c.f. Refs. [23, 24], 



'J a Ap,p) = {p 2 -pI)-(p 2 -m 2 ) , 



(24) 
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that follows directly from Eq. (23). Similarly for the remaining N-2 [a C 2, ....TV— 1] 
dynamically broken O(N) symmetry Nother current matrix elements corresponding 
to 

JU(p>p) = (4> a (p')\Mo)\Mp)) 



(25) 



The identity of the two "gap" , or CDD masses [i\ = M and the concomitant "excess" 
0(2) vacuum symmetry are consequences of the simplicity of the vacuum equation 
(7) that only depends on one 0(N) algebraic invariant [4, 20]. That, in turn, is a 
consequence of the fact that we are dealing with fields in the fundamental irrep. of 
O(N) and the requirement that the Lagrangian (1) be renormalizable, i.e. at most 
of the fourth power in the fields. The assumption of a second v.e.v. (^3) = (tt ) 7^ 
is particularly unrealistic in the Gell-Mann-Levy model (0(N=4) (ft 4 model) [25], 
because of the negative parity of the tt fields: their nonzero v.e.v. would imply 
spontaneous breaking of P and CP "parities". In Sect. 4 we give an example of 
a (f) A model with an internal symmetry that leads to a gap equation with distinct 
mass solutions and potentially exotic NG bound states. 

4 The 't Hooft model 
4.1 Definition of the model 

't Hooft's [26] extension of the linear sigma model Lagrangian reads 



Ah 



tr 



1 2 

- (Ai - A 2 ) [tr (MM*)] - A 2 tr (MM*) 



+ 2k 



where 



e^detM + c.c. 



M = 7y^(£ + *n) 
n = — (tz + tt-t) 



(26) 



(27) 



Eq.(26) is equivalent to the following 



Ah = -[(d„a) 2 + (d^) 2 + (d^rj) 2 + (d^ a ) 
+ 2kcos(9[ct 2 + tt 2 - r/ 2 - a 2 } 
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— 4k sin Q\or] — tt ■ a] 

^1 r 2 , 2 , 2 i 2i2 
o 

- ^[(aa + 7]tv) 2 + (-7T x a) 2 } 



(28) 



which describes the dynamics of the two chiral meson quartets, (a, 7r) and (01,77), 
in this model, and Ai, A 2 , k, 6* are the bare coupling constants. Nonvanishing angle 
6 leads to the explicit (not spontaneous) CP violation in this model, so we set it 
equal to zero. Thus we see that the 't Hooft model consists of two coupled Gell- 
Mann-Levy (GML) linear sigma models [25], one with a light and the other with a 
heavy quartet of mesons. 

Note that the symmetries of various parts of the interaction Lagrangian also 
vary, see App. B: (i) 

\ 1 / 0; A 2 = k = 9 = (29) 

implies 0(8) symmetry, (ii) 

\ 1 / / A 2 ;k = 9 = (30) 
implies 0(4) x 0(2) symmetry, (iii) 

Ai ^0/ A 2 ;k/0 = # (31) 
implies 0(4) symmetry. And the number of NG bosons must change accordingly. 



4.2 The gap equations 

The first set of energy minimization equations (7) reads 

= -v (jiq + 4k) + ^-v v 2 + 3I (m a ) + SI (m a ) + 3lQ(m 7T ) + Io(m v ) 



m„ 



m 



+ 3\2vl (m c 



A 1 



-fiQ - An + yt; v 2 + 3I (m a ) + 3I (m a ) + 3I (m n ) + Io(m v ) 



+ 3A 2 /o ("i c 



-fJ,Q - \K + fy + A 2 ^ V 2 + y [io("V) + 5J (m a ) + 3/ (™7r) + io( m »?)] 



+ A 2 [/ (m (T )+2/ (m w )] 

2 2 -1 1 ^1 2 1 ^1 

™ w = -Mo- 4K +y w +y 

+ A 2 [i"o(^r;) + 2I (m a )] 
2 ,.2 



u 2 + Io(m a ) + 3/ (m a ) + 5/0(771^) + Io(m v ) 



fiQ - 4k + yw 2 + y w 2 + 3io(«v) + 3Io(m a ) + 3/o("i7r) + Io(m v ) 



+ 3A 2 /o("i7r) 



(32) 

(33) 

„)] 
(34) 

(35) 

(36) 



where the divergent integral Jo( m i) is given by Eq. (4). For simplicity in the 
following we use the following short-hand notation a = m^r) = rri^-a = m^;7r = 



10 



m%. Eqs. (32-36) can be solved for 



' - 'I (\r x>-) V " 1,1 ; / 



a-rj (\ 2 -A 2 2 J \ a ~ 7r (J - 7T 



n 7rQ (0) 



J (a) - /o(tt) 



a — 7r 



iW°) 



(a x 2 - 2A X A 2 - 3A 2 2 ) 

a — rj 

1 

(Ai 2 -2AiA 2 -3A 2 2 ) 
whence follows 

two) = 



1 / / 8k \ x t? -8k 
' Ai 1 + A 2 -i 



a — 7T/ a — 7T 




7T — 7/ 

'~k + 7T - 7/ + A 2 (a - rj) U va (0) - 3A 2 (a - ir) 11^(0)) 



n aCT (o) 



Ai (vr - r?) 
io(cr) - 7 (a) 



a — a 



— ^ r f 8k + ^ a - a + A 2 (a - tt) ^(0) - 3A 2 (a - vr) IW0$8) 

Ai (cr - a) \ Ai / 



4.3 The Nambu-Goldstone theorem in the isotensor pseu- 
doscalar sector 

Note that there are 8 x 8 = 64 possible initial or final two-body states here. Thus 
there are 64 x 64 = 4096 possible channels, but only 4 x 7 = 28 distinct pairs 
of particles, or equivalently at most 28 possible NG bosons. The last statement 
holds under the proviso of exact 0(8) symmetry being broken to a discrete (non 
Lie) symmetry, however. Otherwise there are fewer than 28 NG bosons, and when 
0(8) is explicitly broken down to 0(4) x 0(2), or 0(4) by one of the terms in the 
Lagrangian Eq. (26), there are even fewer than that. This fact tells us that many 
channels must be coupled and that the residual symmetry plays a crucial role in 
this coupling. We shan't look at every possible channel in this paper, but rather 
concentrate only on the exotic ones, in this case the isotensor, so as to show the 
existence of pure bound state NG bosons without CDD admixtures. 

We may use the residual vacuum symmetry, e.g. the 0(3) isospin invariance to 
split this 64 x 64 matrix equation into six invariant subspaces: three flavour channels 
[(a) isoscalar; (b) isovector, and (c) isotensor] of either parity. In the two-body, or 
Bethe-Salpeter (BS) equation for the four-point Green functions Dij(s), the indices 
i,j denote the isospin of the two-body initial and final states, respectively. 
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The negative parity (pseudoscalar) isotensor two-body equation is a single- 
channel one and straightforward to solve, see Eq. (12). We look at the zero CM 
energy P = function y 7rQ (0)Il 7ra (0); we use 

Vy=V = A 1 + A 2 (39) 

and Eq. (37) to obtain the final result 

Now set A2 = k — > and find 

, lim V&=*> (0)11™ (0) = 1 • (41) 
The propagator Eq. (12) evaluated at zero momentum can be written as 

dV=2)(o) = 9 ™ = v ™(°) (42) 

™ [ ' -ml a 1-^(0)11^(0)' 1 > 

whence it follows that 



m ira 



1 ' Ai ( 1 I + A 2 



Ai + A 2 VAi - 2AiA 2 - 3A 2 / V V a — it J a - tt 
a — it A - 



v 

= 0(X 2 ) + 0(k) (43) 

Thus we see that the effective (pseudo) NG boson mass in this channel is propor- 
tional to A 2 , and/or k, the two 0(8) symmetry breaking parameters, q.e.d. 



5 Summary and Conclusions 

In summary, we have: 1) proven the NG theorem in the variational Gaussian wave 
functional approximation to the O(N) symmetric (f) 4 model when the symmetry 
of the ground state is a proper subgroup of O(N-l); 2) proven conservation of 
Nother currents corresponding to the dynamically broken symmetries; 3) proven 
the same NG theorem in the exotic isotensor channel of the 't Hooft model in the 
limit A 2 = k — > 0. The NG bosons are massless bound states of two massive 
constituents. We emphasize that our proofs do not depend on the specific values of 
the bare parameters, or of the cutoff in the theory, so long as the system is in the 
spontaneously broken phase with appropriate symmetry. 

We should like to put these results into their proper logical and chronological 
setting. The variational method in quantum field theory (QFT) is based on the 
Schrodinger representation and goes by the name of Gaussian approximation to the 
ground state wave functional. This method, in its various guises, was pioneered 
by Schiff, Rosen and Kuti [6, 7, 8] in the 1960's and 70's, and later revived and 
elaborated in the 1980's by Barnes and Ghandour [9, 13], and by Symanzik [10] and 
by Consoli, Stevenson and collaborators [11, 12, 16]. Related formalisms based on 
effective potentials and other functional methods were discussed in Ref. [8, 16, 12] 
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and references cited therein. Most of these studies addressed the <p A scalar field 
theory that is also the prime example of the Nambu-Goldstone (NG) theorem [1, 2, 
3, 4, 5], an exact result in the O(N) symmetric <fi A scalar theory with spontaneous 
internal symmetry breaking. 

The NG theorem was first shown not to be satisfied by the solutions to the mass, 
or "gap" equations in the Gaussian approximation by Kamefuchi and Umezawa in 
1964 [14], practically simultaneously with the general proofs of the NG theorem 
[3, 4, 5]. This fact was subsequently rediscovered several times [15, 16] and this 
unsatisfactory situation persisted until 1994 [17]. Various conjectures as to the rea- 
sons for this failure and as to potential remedies were advanced during this period 
of time. It was first shown in Ref. [17] that this apparent breakdown of the NG 
theorem is not an intrinsic shortcoming of the Gaussian wave functional approxima- 
tion, but rather a consequence of incomplete previous analyses. In other words, the 
NG theorem is satisfied, but the NG bosons are not excitations of the elementary 
scalar fields, as initially expected. Rather, they are massless bound states of two 
massive elementary scalar excitations, in close analogy with Nambu's [1] proof of 
the NG theorem in a self-interacting fermion theory. NG bosons are solutions to the 
two-body (or Bethe-Salpeter [BS]) equation in the Gaussian approximation. That 
equation was only rarely considered in the literature [7, 9], and never before Ref. 
[17] in the context of spontaneous symmetry breaking of purely bosonic models. 

In this light the result is simple enough to understand, yet it drew strong, albeit 
unpublished criticism and affirmation [27]. Perhaps the underlying reason for the 
misunderstanding by some was the implication of the proof that the (J) 4 scalar field 
theory could have bound states, which, as "everybody knew" [27], disagrees with 
various "rigorous no-go" and "triviality" theorems in the same theory [28] 1 . The 
said theorems hold only in the limit of an infinite cutoff, however, in which the 
Gaussian approximation also becomes trivial [11]. For finite cutoffs, on the other 
hand, this is a nontrivial theory that may contain bound states. 

Soon after the first proof in Ref. [17] it was also shown along the lines of Ref. 
[3] that the NG theorem also follows from the Gaussian effective potential, Ref. 
[30], but that proof did not shed much light on the mechanisms that made the 
NG bosons come about. Only later, in Ref. [31, 32], it was explicitly shown how 
this formal proof relates to the Gaussian two-body equations of motion. Another 
source of confusion was the apparent doubling of degrees of freedom, at least in some 
"flavour" (internal symmetry) channels, viz. the existence of massive "elementary" 
and massless bound states in the same channel. This problem was resolved in Ref. 
[31], wherein the Kallen-Lehmann spectral function was calculated in appropriate 
channels of the model. This spectral function clearly shows the presence of massless 
NG states and the absence of the massive single-particle excitations. That also 
constitutes a proof of the NG theorem along Gilbert's lines, Ref. [5], within the 
Gaussian wave functional approximation. Thus we have confirmed all the well 
known proofs of the NG theorem in the Gaussian approximation. 

The NG theorem is the simplest example of a Ward-Takahashi identity, which 
follows from the underlying internal symmetry of the (j) 4 model. Ward-Takahashi 
identities typically relate (n-l)-point Green functions to n-point functions and/or 

7 For more recent results comparing constructive QFT to the Gaussian approximation in 1+1 dimen- 
sional field theories, see Ref. [29] 
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matrix elements of Nother currents. These identities were developed by Lee in the 
linear sigma model at the perturbative one loop level [23], and by Symanzik for 
arbitrary orders of perturbation theory [24], so we shall call them the Lee-Symanzik 
[LS] identities. 

The exact, i.e. nonperturbative Green functions satisfy an infinite set of coupled 
integro-differential equations called the Schwinger-Dyson equations [19]. The itera- 
tive/perturbative solutions to the SD equations form (infinite/finite) sets of Feyn- 
man diagrams. If one decouples the SD equations for higher order Green functions 
from the lower order ones (in popular jargon, if one truncates the SD equations), one 
may obtain tractable equations and find their solutions that sum infinite, albeit in- 
complete sets of Feynman diagrams. It has been known at least since 1980, Ref. [9], 
that the Gaussian approximation to the unbroken symmetry </> 4 theory corresponds 
to one such truncation of the SD equations. But, truncated SD equations need not 
obey the conservation laws of the original SD equations of motion, i.e. LS identities 
may be violated by the truncation. To our knowledge, no proof of LS identities had 
been given for truncated SD equations, i.e. for infinite classes of diagrams in the 
bosonic linear sigma model prior to Ref. [17], although a similar proof had been 
given by Nambu and Jona-Lasinio in their fermionic model some 30 years before 
[1]. Thus we have shown that the Gaussian functional approximation constitutes a 
closed, self-consistent symmetry-preserving approximation to the Schwinger-Dyson 
equations. 
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A The gap equations with two v.e.v.s 

The first set of energy minimization equations (7) read 



Ar 



v 2 + 3hI (M) + hl (m) + (N- 2)hh(n 2 ) 



A = A [v 2 + hI (M) + 3M ( Pl ) + (N - 2)M (p 2 ) 

v 2 = (M 2 + (<Pif = <0> 2 
<&> = i = 2,...,JV-l , 



(44) 
(45) 

(46) 
(47) 



where the divergent integral Io(rrii), Eq. (4) is understood to be regularized via an 
UV momentum cut-off A, either three-, or four dimensional. The second set of gap 
equations, (7) read: 

M 2 = -ii 2 + A [2{4> ) 2 + (0) 2 + 3hI (M) + ni (m) + (N- 2)M (/x 2 )] (48) 
li\ = - f j 2 + \ \2{<p 1 ) 2 + {cf>) 2 + hI (M) + 3hIo(iii) + (N-2)hIo(v2)}m 



+ A (<£> 2 + ftio(M) + hlobn) + iV?i/o(^2 



(50) 
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Upon inserting Eqs. (44, b) into Eqs. (48, b), the following coupled "gap" equations 
emerge: 

M 2 = 2A o (0o) 2 + 2A o ?i[/o(M)-/ o (/i 1 )] = 2A o (0o) 2 (51) 

n\ = 2A (<M 2 -2AoM/ (M)-/ (/i 1 )] = 2A o (0i) 2 (52) 

M 2 -fij = 2Ao((</>o) 2 -(0i) 2 + 2A o + ?i[/o(M)-/o( m )]) (53) 

$ = 2X h [/ (/i 2 ) - I (M)} = 2X h [/o(M2) - /o(Mi)] • (54) 
Note that these equations lead to 

/ (M)-/ (/ii) = (55) 

M 2 -^ = 2A ((0 O ) 2 - (^i) 2 ) = , (56) 

which, in turn have at least one solution: \i\ = M and (4>o) 2 = (0i) 2 5 that are 
solutions to a single nontrivial gap equation 

M 2 = 2A (<M 2 (57) 
& = 2\ h[I (p 2 )-I (M)] , (58) 

with two unknowns one of which is kept fixed, c.f. Ref. [18]. This gap equation 
has been solved numerically in Ref. [18]: it admits only massive solutions M > 
H > 0, however, for real, positive values of Xo,l^o an d real ultraviolet cut-off A in 
the momentum integrals Io(rrii). In other words the "would be NG boson fields" 
(01,2,.. .,w-i) excitations are all massive (fj, > 0) in the MFA. This looks like a 
breakdown of the NG theorem in this approximation, but, as discussed in Ref. 
[17, 18], there is a solution by way of the two-body (Bethe-Salpeter) equation. 

B Symmetries of the 't Hooft model 

The two field quartets, (ct,7t) and (at, if), have different "crural" 0(4) = 0(3) x 
0(3) ~ SU L (2) x SU R (2) ' 

<5 5 cr = (3 ■ 7T (59) 

5 5 r] = -/3ct (60) 

5 5 a = Pr) (61) 

5 5 7T = -Pa , (62) 

isospin 

5a = (63) 

5r] = (64) 

5a = —e x ct (65) 

(57T = — £ X 7T , (66) 
and £7a(1) — 0(2) transformation properties 

5° 5 a = fa (67) 

5lv = -Per (68) 

5% a = (3tv (69) 

^g-n- = -Pa . (70) 
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The Lie algebra 0(4) has two Casimir operators, but there is only one invariant 
in the (5,5) representation with one meson quartet, viz. (i) [a 2 + 7r 2 ], whereas 
with two quartets one has three invariants: (i) above; (ii) [rf + a 2 ], and (hi) [r]a — 
a ■ 7r]. Any odd power of the third invariant (hi) violates CP. Even without the 
third invariant the algebraic structure of the Lagrangian is rich enough to allow 
for multiple vacua [solutions to the energy minimization equations (6)] even in 
the (first) Born approximation. For example, equations (6) applied directly to 
the 'tHooft interaction (26) in the Born approximation allow two nonzero v.e.v.s 
[(f)oB = vqb ¥= and (a 3 ) B = «ib ¥= 0]: 



^OB + «1B + X 2 V 1B 



9 Al 

-f4 -4k+ — 




0. 



Their solutions are 



(o")ob = "Wob = 

(«3)0B = Vib = 



Ai + A2 A2 
Ai + A2 A2 



(71) 
(72) 

(73) 
(74) 



leading to a nontrivially broken ground state. Note, however, that (i) A2 — > is 
a singular limit point, i.e. these vacua are not continuously connected with the 
more conventional vacua at A2 = 0, and (ii) the two kinds of vacua/energy minima 
coincide whenever k = and A2 7^ 0. Hence it is possible for the system to be 
in an unconventional vacuum even with infinitesimally small A2 and k when they 
vanish with a fixed nonzero ratio. This makes it plausible that unusual vacua may 
also appear in this model when the symmetry is broken dynamically, i.e. by "loop" 
effects. This doubling of vacua is a consequence of multiple (two) independent 
algebraic invariants in the Lagrangian (26) which in turn led to multiple (two) 
v.e.v.s, in agreement with the general theory [4, 20]. 
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